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Abstract. We simulate short positive and negative streamers in air at standard 

temperature and pressure. They evolve in homogeneous electric fields or emerge 

from needle electrodes with voltages of 10 to 20 kV. The streamer velocity at given 

streamer length depends only weakly on the initial ionization seed, except in the case 

t/^ . of negative streamers in homogeneous fields. We characterize the streamers by length, 

,-H I head radius, head charge and field enhancement. We show that the velocity of positive 

•"T^' streamers is mainly determined by their radius and in quantitative agreement with 

Y^ , recent experimental results both for radius and velocity. The velocity of negative 

• i-H ' streamers is dominated by electron drift in the enhanced field; in the low local fields 

C/2 ' ...... ... 

^.- . of the present simulations, it is little influenced by photo-ionization. Though negative 

,^ I streamer fronts always move at least with the electron drift velocity in the local field, 

, '"*' this drift motion broadens the streamer head, decreases the field enhancement and 

ultimately leads to slower propagation or even extinction of the negative streamer. 
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1. Introduction 

Atmospheric pressure corona discharges are widely used in technology. Streamers, which 
are the basic building blocks of these discharges, focus a large part of the energy of the 
reactor into a small volume. As positive streamers emerge from pointed electrodes at 
lower voltages than negative ones [H E], recent investigations have largely focused on 
positive streamers. However, negative streamers are clearly present in many natural 
phenomena of atmospheric electricity like lightning, sprite discharges etc. P, [3], HI [5]. 
Furthermore, modern high voltage supplies easily create negative streamers [2l[6l[7], and 
they are very promising for disinfection applications [8] if electrical matching problems 
can be overcome. An experimental study of positive and negative streamers in air at 
standard temperature and pressure in a wide voltage range is available in [2]. 

The simulation of positive streamers in three spatial dimensions with cylindrical 
symmetry meanwhile is based on a large body of research. Pioneering work was done 
by Wang and Kunhardt [9] and by Dhali and Williams [lO] . The use of more complete 
and realistic plasma-chemical models [11], better modeling of the electrode geometry 
[T2I [T3] and an efficient calculation of the non-local photo-ionization source [Ml [15] 
have finally allowed simulation and experiment results to converge within a narrow 
range. Pancheshnyi et al. [16] were able to predict the mean streamer velocity at varying 
pressures within a range of around 25% and thus question the role of photo-ionization 
versus fast electron detachment in repetitive positive streamer discharges in air |17j . 
Also remarkable was the reproduction of experimental results [18] of streamers in long 
gaps of 13 cm performed in |19j . 

In early work, the non-local ionization mechanism through photons was replaced 
by background ionization, and positive and negative streamers looked fairly similar. An 
example of a simulation of a double-headed streamer that is completely dominated by 
the assumption of the initial ionization distribution, can be found in [20] . Since photo- 
ionization was introduced as a non-local ionization mechanism in air to explain the 
propagation of positive streamers, precisely three groups of authors have investigated 
negative streamers in air with photo-ionization. Babaeva and Naidis have compared 
positive and negative streamers emerging from pointed electrodes in a short paper in 
1997 [21], Liu and Pasko have investigated doubled-headed streamers in homogeneous 
fields [22| [23] , and the present authors have studied the infiuence of photo-ionization on 
propagation [15] and interaction [21] of streamers of both polarity. 

The present paper is devoted to a systematic study, characterization and 
comparison of positive and negative streamers in ambient air. It is organized as 
follows. In section [2] we describe our model. Section [3] treats double-headed streamers 
in homogeneous fields, their dependence on the ionization seed chosen as a starting 
point for the simulations and their basic mode of propagation. Section [H treats positive 
and negative streamers emerging form needle electrodes; streamers of both polarity are 
characterized by velocity, field enhancement, head radius and head charge; characteristic 
differences are found. Their velocities are dominated either by the head radius for 
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positive streamers or by the enhanced field for negative streamers. Section \5\ shows a 
convincing comparison with the experiments in [2]. Finally, we summarize our main 



results in section El Appendix A contains the Charge Simulation Technique for the 



needle electrode, and Appendix B the definition of the quantities used for the electro- 



dynamic characterization of the streamer. 

2. Model 

We use a fiuid model of air that contains electrons and six species of ions: N^, O^, 
N^, Oj, 0^N2 and O^. While electrons diffuse and drift in a self-consistent electric 
field, the ions due to their much larger mass can be approximated as immobile. We 
consider 15 reactions among the species, taken from [16] plus photo-ionization which 
we approximate as described in [15]. All our simulations are performed at standard 
pressure and temperature. We assume cylindrical symmetry of the streamer and we 
solve the model numerically by means of adaptively refined grids [25] . 

We performed numerical experiments in two geometries: the simplest one is defined 
by plane parallel electrodes with fixed electrical potential, therefore we impose Dirichlet 
boundary conditions to the potential. A more complex geometry is given by a needle 
electrode, which we simulate by means of a simplified version of the Charge Simulation 
Technique, as detailed in [Appendix A 



To facilitate the physical interpretation of the simulation results we provide, 
besides some direct results, a number of derived quantities to characterize the streamer 
evolution. The precise definition of these quantities (length L, velocity v, enhanced field 
Emax, radius R and charge in the streamer head Q) is provided in Appendix B The 



choice of these particular quantities is guided by the aim to develop an electro-dynamic 
model of streamer heads and streamer channels in the future; qualitative models of this 
type have been suggested in the past Russian literature [261 1211 [281 E]. However, these 
models require further improvements as was discussed in [29] . 

3. Double-headed streamers in homogeneous background fields 

3.1. Motivation 

We first study streamers in homogeneous background fields. Most streamer experiments 
are performed in needle-plane [30l [STl [2], wire-plane or wire-cylinder geometries [7], 
some also between planar electrodes with a protrusion [T8l [32] : here the protrusion is 
used as the inception point of the discharge. Experiments of streamers between planar 
electrodes created with a laser were difficult to interpret [33]. The last two experiments 
were particularly designed to study streamers in homogeneous background fields while 
the inception form curved electrodes is much easier. But even if one electrode is a 
needle or a wire and the other one a plane, once the streamer tip approaches the planar 
electrode, the background field is again well approximated by a homogeneous field. 
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Figure 1. Electric field of a double-headed streamer extending between two planar 
electrodes, plotted at equal time steps 1.2 ns. The negative front is propagating 
upwards, the positive moves downwards. The number of particles in the initial seed 
is 6 • 10^". Note that the lateral borders of the figure do not correspond to the 
computational domain. 



Finally, the electric field responsible for sprite discharges, located between a charged 
thundercloud and the ionosphere, is rather homogeneous. 



3.2. Simulations and dependence on the initial ionization seed 

A streamer discharge in a homogeneous field is initiated by a localized ionization seed 
created, e.g., by a cosmic particle shower or by an electric field inhomogeneity around a 
suspended particle. If the field is above ~30 kV/cm, a small seed first can undergo an 
avalanche phase where the ionization level increases exponentially and then, once there 
are significant space charge effects, it reaches the streamer regime [M]. If the field is 
below threshold, a seed of finite size is required that after drift separation of charges 
immediately enters the streamer phase. 

To investigate whether the initiation mechanism in a field above threshold has a 
lasting effect on the propagation of a streamer, we run several simulations with different 
initial seeds. The background field was Eq = 50kV/cm and the seeds had a spherical 
Gaussian profile with an e-folding radius of 74 /^m located at the center of a 14mm-gap, 
creating a double-headed streamer. The temporal evolution of the streamer shape is 
illustrated in Fig. [1] for a particular seed. Panel a of Fig. [2] shows the position z of the 
two ionization fronts as a function of time t; it shows that both fronts accelerate in time, 
and that the fronts with the particle rich seeds are ahead in evolution to those that start 
with a weaker seed; they propagate faster at any particular time both on the positive 
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Figure 2. a) Position z{t) of the negative (upper) and positive (lower) fronts of a 
double-headed streamer as a function of time t. The streamer evolves between two 
planar electrodes with a background field Eq = 50 kV/cm as shown in Fig. [TJ The 
distance between the electrodes is 14 mm. The different colors correspond to different 
numbers of initial electrons and NJ ions ranging from 6 • 10^ to 6 • 10^° as indicated 
in the figure (the seed is initially neutral), b) The same data, now plotted as the 
absolute values of velocities v{z) of the negative (right) and positive (left) fronts as a 
function of the front location z. Note that each front propagates in a different direction, 
but the absolute value of the velocity is plotted to help the comparison. The dotted 
horizontal line indicates the velocity v*{Eq) of a planar negative ionization front in the 
background field Eq. 



and on the negative side. 

The front velocity v as a function of time t is simply the derivative of the curves in 
panel a, and v(t) strongly depends on the size of the initial seed both for the positive 
and for the negative front. On the other hand, panel b of Fig. [2] shows the absolute value 
of the velocity f as a function of position z — this is the observable typically measured 
in the experiments. Here the curves v{z) for the positive fronts essentially overlap for 
different seeds while those for the negative fronts don't: the negative streamers with the 
largest seed have the largest velocities v{z) as a function of head position. 



3.3. Discussion of inception and propagation of positive and negative streamers 

We now discuss the physical mechanisms causing the different dependencies of v{z) on 
the initial seed. The horizontal dotted line in panel b of Fig. |2] indicates the velocity 
v*{Eo) where 

v*{E) = \E\+2^/ir\E\~MlE\) (1) 

is the velocity of a planar negative ionization front in a field E in the absence of photo- 
ionization in dimensionless units [361 ES]; it is given by the local electron drift velocity 
\E\ augmented by the combined effect of electron diffusion D and impact ion ionization 
q;(|£^|). The velocity v*{Eq) evaluated in the background field Eq is a lower bound for the 
velocity of a negative ionization front with field enhancement and with photo-ionization 
as indeed can be seen in the figure. Obviously, the negative front will always propagate 
at least with the electron drift velocity in the local electric field. If the initial seed is 
stronger, field enhancement will build up faster while the front is already in motion, and 
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at each position z, the front is faster than for a weaker seed. The positive front, on the 
other hand, has no lower bound for its velocity. The positive discharge side stays at rest 
until photo-ionization has built up a sufficient electron concentration for the streamer 
to start. This means that the inception time now strongly depends on the seed, but 
once the positive streamer propagates, it does it with a similar velocity v as function of 
the position, rather independently of the seed. 

In view of the experimental results, a most interesting question is the comparison 
of the velocities of positive and negative streamers. If we fix the time elapsed after the 
seed of the double-headed streamer is created, the negative streamer is faster than the 
positive one. If, on the other hand, one compares the velocities at a fixed distance from 
the initial seed, the picture is different: for small distances, the negative streamers are 
faster, but they are overcome by the positive ones at larger distances from the initial 
seed. Note also that the differences between positive and negative streamers become 
the smaller the larger the initial number of particles. 

The only other studies of double headed streamers with photo-ionization in a 
homogeneous electric field are performed by Liu and Pasko [221 123]; in particular, table 2 
in [22] shows that the characteristics of streamer propagation, namely field enhancement 
in the streamer head and field screening in the streamer interior are stronger on the 
positive streamer head, and that the positive streamer is faster. This is found at air 
models applicable to a height of 0, 30 and 70 km in the atmosphere. 

The larger velocity of the positive streamers is surprising if one takes into account 
that for identical field enhancement and identical electron distribution in the leading 
edge of the ionization front, the negative front will always be faster [36] as electron drift 
supports propagation of negative fronts and acts against it for positive fronts. However, 
inspection of Fig. [1] shows that the positive streamer is more focused and the field at its 
head is more enhanced. Ultimately, electron drift leads to a "dilution" of head focusing 
and field enhancement in negative streamers and makes them run slower at a given 
distance form the ionization seed. 

4. Streamers in needle-plane geometries 

The study of streamers disconnected from electrodes in homogeneous fields already 
gives a good qualitative insight into many of their properties. However, laboratory 
experiments and engineering applications are mostly done in inhomogeneous fields and 
streamers emerge from pointed electrodes. We therefore here study positive and negative 
streamers emerging from a needle electrode and propagating towards a planar electrode. 
This electrode configuration was implemented by means of the Charge Simulation 



Technique as described in Appendix A 
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Figure 3. Front propagation of negative and positive streamers for different seeds 
between needle-plane electrodes. The gap is L = 11.5 mm long with an applied voltage 
of 23 kV; the needle parameters are Lneedie ~ 2.3 mm, Rneedie = 0.26 mm. Different 
colours correspond to different number of particles in the initial seed. 



4-1. Weak dependence on the initial conditions 

In order to study how the streamer behavior depends on parameters, we performed a 
number of different simulations. The first observation is that in contrast to the case of 
homogeneous fields, in inhomogeneous fields the initial seed affects the propagation of 
the streamer only slightly, even when considering the streamer velocity t" as a function 
of time t. This is illustrated in Fig. [3] where the positions z of streamers as a function 
of time t are shown for several initial seeds — this plot corresponds to panel a in Fig. O 
The seeds are placed at the top of the gap and have spherical Gaussian profiles with 
radius 92/_im and between 6 ■ 10^ and 3 ■ 10^ particles. 

For positive streamers in needle-plane geometries, this effect was found before by 
Pancheshnyi et al. [37]. The reason is probably that in the high field region near the 
needle electrode, the seed grows very rapidly compared to the slower evolution in the 
lower fields away from the needle. 



4-2. Simulations of positive streamers 

We now analyze positive streamers in the needle-plane geometry in more detail. We 
run simulations of positive streamers in a short gap of 7 mm gap at a voltage of 
10.5 kV, 14 kV and 21 kV. The needle electrode had a radius Rneedie = 0.2 mm and 
a length Lneedie = 2 mm. An example of the evolution of the streamers is shown in 
panel a of Fig. HI An array plotting the relationship between each pair of the derived 



quantities defined in Appendix B for each of our runs is provided in Fig. [51 The purpose 
of this array representation is to find consistent relationships between the streamer 
characteristics. The overall evolution shown in Fig. [5] can be summarized as follows: 
as the streamer advances, it becomes thicker and faster and the total charge in the 
streamer head increases while the enhanced electric field decreases. 
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a) 





Figure 4. Electric field created by a propagating streamer in a needle-plane geometry 
plotted at equal time steps. Note that the computational domain extends beyond the 
lateral borders of the plot, a) A positive streamer at time steps of 2.7 ns in a voltage 
of 14 kV, b) a negative streamer at time steps of 4.5 ns also in a voltage of 14 kV. 



4-3. The velocity of positive streamers 

While for a negative streamer, the velocity v cannot become smaller than the electron 
drift velocity I-Emaxl in the locally enhanced field, for the positive streamer the 
velocity v increases while the field enhancement I-Emaxl decreases. In the experimental 
investigation [2], the completely empirical relation 

0.5 d^ 

V ^ 2 

mm ns 

between velocity v and diameter d of positive streamers in air at standard temperature 

and pressure was fitted to the experimental data (in Fig. 6b of [2]). In Fig. [5], the 

plots of f as a function of radius R also lie more or less on one line. In Fig. [6] we 

therefore compare the v{R) plot of our simulation results with the empirical fit ([2]) 

to the experiments. The plot does not contain any fitting parameters and shows a 

close resemblance between experiment and simulation. When interpreting the figure, it 

should be noted that the simulations measure the geometrical radius of the space-charge 

layer, also called "electro-dynamic radius", while an experiment measures the visible, or 

radiative, radius. There can be a significant difference between both measures; in [16] 

it is estimated that the electro-dynamic radius is about twice the radiative radius. 

A positive streamer propagates due to the photo-ionization in front of its head. 

Comparing the head radius to the photo-ionization lengths [2l], the hypothesis of 

Kulikovskii [38] that the streamer radius would be determined by the photo-ionization 

length can be clearly discarded. The photo-ionization absorption is fitted by two 

lengths scales, one of them is essentially negligibly small, the other one much larger 

than the head radius [IHl El] . Free electrons created by photo-ionization are therefore 

available throughout the region where the electric field is above the ionization threshold. 

Therefore the streamer velocity will be mainly determined by the size of the region where 

further electron multiplication by impact is efficient. This region, in turn, is roughly 
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Figure 5. Characterization of positive streamers by their main (electro-)dynamical 
variables, namely time t, length L, maximal field Emaxi velocity v, head radius R 
and charge in the streamer head Q; these quantities are defined in [Appendix B[ The 
picture shows each variable as a function of each other variable; the plot with variable 
A in the X axis and variable B in the Y axis is located in the column and row which 
have, respectively, A and B in their diagonal cells. Hence, every figure appears twice 
with its axes swapped. Red circles correspond to a voltage of 10.5 kV, black crosses 
to 14 kV and gray triangles to 21 kV. Note that all plots are qualitatively similar for 
different applied voltages. Remarkably, the plots of streamer velocity v as a function 
of radius R nearly fall on the same line for different voltages. 
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Figure 6. Relationship between velocity v and radius R of positive streamers. The 
symbols are simulation results as presented in Fig. [51 the continuous line represents 
the empirical fit ^. 



determined by two factors: the enhanced electric field and the electro-dynamic radius 
of the streamer. In Fig. [5] we see that the radius varies much more than the enhanced 
field, which explains the qualitative relation between velocity and radius. 

4.4- Simulations of negative streamers 

We investigated negative streamers also in a 7mm-gap with applied voltages of 10.5 kV 
and 14 kV. The electrode geometry was the same as for positive streamers. The 
evolution of the spatial structure is illustrated in panel b of Fig. HI As in the case 
of the homogeneous field in Fig. [H the negative streamer is broader and the field is less 
enhanced. 

Figure [7] shows the relationship between relevant (electro-) dynamic quantities 
during the evolution of these two streamer simulations. The streamer in the lower 
voltage (red symbols) becomes slower, the enhanced field eventually drops below the 
threshold for impact ionization, the head charge disappears and the radius diverges. 
At that instant of time, the remaining electrons from the streamer head continue to 
drift towards the planar electrode, but the impact ionization ceases to be efficient 
and the streamer mode of propagation stops. This is probably the generic way how 
a negative streamer extinguishes, quite different form the one reported for positive 
streamers in [39]. The streamer in the higher voltage undergoes a similar intermediate 
evolution. However, eventually the proximity of the planar electrode again enhances the 
field and the streamer reaches the electrode. 



4-5. Velocity of negative streamers 



As already discussed above and in section [331 negative streamers propagate not only due 
to photo-ionization but also due to electron drift. For the velocity, this implies a stronger 
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Figure 7. Characterization of negative streamers by their main (electro-)dynamical 
variables, the presentation is the same as in Fig. [5] for positive streamers. All 
simulations were performed in the same geometry as for the positive streamers, with 
applied voltages 10.5 kV (red circles) and 14 kV (black crosses). Note that in the low 
voltage case, the enhanced field ceases to be strong enough to sustain the streamer 
propagation. 
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Figure 8. Relationship between velocity v and enhanced field Emax for negative 
streamers. The symbols arc simulation results as presented in Fig. [71 the continuous 
line represents the velocity of a planar front without photo- ionization v*{Ejnax) ©• 



dependence on the enhanced field, that thus overcomes the dependence on the radius. In 
fact, comparison of Figs. [5] and [7] shows that the velocity of positive streamers increases 
with radius while the field enhancement decreases; the velocity of negative streamers, 
on the other hand, increases with field enhancement while the radius decreases. Indeed, 
Fig. [7] shows a very clear correlation between velocity v and field enhancement Emax for 
negative streamers that we now analyze further. 

Actually, one can compare the actual velocities with the velocity v*{Emax) where 
v*{E) is given in Eq. ([1]). The velocity v*{Emax) is the velocity of a planar fully relaxed 
negative streamer front in a field E^ax when the effect of photo-ionization is neglected. 
Fig. M shows the simulation data for v{Emax) from Fig. [7| and the function v*{Emax) 
for comparison. The coincidence is strong. Deviations mainly come from the fact 
that the analytical equation is a lower bound to the actual velocity as photo-ionization 
is neglected (for a comparison of simulation data without photo-ionization with the 
analytical formula we refer to [iQl ST]). 

One should remark here that the background electric field under which the streamer 
propagates most of the time is quite low and hence one does not observe a buildup of 
the ionization level in front of the streamer as in the high-field case of [I5]. Therefore, 
we do not observe a transition to a regime dominated by photo-ionization. 



4.6. Comparison of positive and negative streamer simulations 

We here directly compare the propagation of a positive and a negative streamer in a 
longer gap of 11.5 mm than considered above in needle-plane geometry. For both positive 
and negative streamers the applied voltage is l^ = 23 kV. The radius of the needle is 
Rneedie = 0.26 mm and its length is Lneedie = 2.3 mm. The position as a function of time 
was already presented in Fig. [3] where it was shown that the evolution depends only 
very weakly on the initial ionization seed. The figure shows that the negative streamer 
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Figure 9. Electric field on the streamer axis at equal time steps of 3 ns for negative 
(upper) and positive (lower) streamers propagating to the left. The number of particles 
in the initial seed is 6 ■ 10^. 



initially is faster, but it is soon overtaken by the positive one. 

In Fig. [9] we show the spatial profiles of the electric field on the streamer axis for a 
number of time steps; the streamers are propagating to the left. The electric field at the 
positive streamer heads is much more enhanced than on the negative ones. This larger 
field enhancement is due to the smaller radii of the positive streamers that consecutively 
propagate much faster. Also the field inside the streamer channel is screened less for 
negative streamers. We note that in the only other comparable simulation of positive 
and negative streamers by Babaeva and Naidis [21], the field inside the negative streamer 
channel is higher as well, but Fig. 6 of that paper shows that their negative streamers 
are faster than their positive ones, though there is also a consistency problem between 
their Figs. 6 and 7, and positive and negative streamer velocities are not compared in 
the text. 



5. Comparison with experiments 

Experiments [71 [2] show that positive and negative streamers in ambient air driven by 
voltages above 60 kV behave qualitatively similar. On the other hand, there are major 
differences below 40 kV [2]. Positive streamers form at lower applied voltages, they 
are faster, longer and thinner. Our simulations at voltages between 10.5 and 21 kV in 
shorter gaps reproduce all these features. 



5.1. Inception 

While the full inception process in interaction with the electrode needle surface is not 
part of the present simulations, we observe that positive streamer inception is not very 
sensitive to the initial ionization seed while the negative streamer formation depends on 
it, at least in homogeneous fields. 
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5.2. Velocity 

Close to the needle electrode, the electrons of the negative ionization seed drift outwards 
in the local field and are faster. However, just the lack of outward drift motion in the 
positive seed leads eventually to a larger field enhancement and ultimately to a faster 
propagation of the positive streamer at the same distance from the electrode. 

Experimental measurements of the velocity of positive streamers are very well fitted 
by the empirical equation ([21) that relates their velocity to their radius. Figure [6] shows 
that this equation also fits our simulation results quite well without any fit parameter. 

5.3. Diameter 

The negative streamers are thicker and less focused, both in simulations and in 
experiment. The minimal diameter of positive streamers in our simulations is about 
0.2 mm, identical to the minimal diameter reported in experiments [301I3T1I2]. It should 
be noted, however, that the definition of the radius might differ between experiments 
and simulations as discussed in section 14.31 

5.4. Length and extinction 

In a potential of 10.5 kV, the negative streamers extinguish after less than 2 mm, 
while the positive ones reach the planar electrode at 7 mm distance. In fact, in the 
experimental paper [2], discharges of 2 mm length are not called streamers, and the 
extinction of these very short negative discharges is in agreement with experiment. At 
14 kV, our simulated negative streamers do reach the planar electrode, but they are 
helped by a strong initial ionization seed and a short gap. Simulations in longer gaps 
are in progress. 

6. Conclusion and outlook 

We have studied the propagation of double-headed streamers in a homogeneous field 
and the inception and propagation of positive and negative streamers emerging from 
needle electrodes. We have shown that for spatially concentrated ionization seeds 
containing from about 10^ to about 10^'^ electron-ion pairs, the streamer velocity at 
a given streamer length depends only weakly on the seed, except for the case of negative 
streamers in homogeneous fields. We have found qualitative and quantitative agreement 
with experiment as summarized in section [51 

We have shown that the relations between velocity v, radius R, field enhancement 
Emax and head charge Q that characterize the propagation of a streamer, differ 
qualitatively between streamers of different polarity. The velocity of a positive streamer 
in air is mainly determined by its radius, in accordance with the empirical fit formula ([2|), 
while that of a negative one is dominated by the enhanced electric field and well 
approximated by v*{Emax) where v*{E) is the velocity of a planar negative ionization 
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front in the absence of photo-ionization. These statements on the simulations hold in 
the voltage range of 10 to 20 kV in gaps of 7 or 11.5 mm length. Longer gaps and higher 
potentials will be investigated in the future. 

The characterization of the streamer head by velocity, radius, field and charge is a 
first step towards an electro-dynamical characterization of streamer head and channel. 
Such models were already sketched in [261 1211 [Ml [3], however, they need to take care of 
the polarity dependence of streamers, and they have to be made consistent with charge 
conservation |29j . 

The fact that in experiments [2] positive streamers moved faster than negative 
ones, initially was quite puzzling from a theoretical point of view. As electron drift acts 
against positive and in favor of negative streamers, a streamer with identically formed 
space charge layer in its head will always move faster, if it is negative. However, the 
simulations show that the electron drift does not necessarily help the negative streamer 
to propagate; rather the outwards drift motion that is essentially linear in the field, leads 
to a growth of the head radius and a subsequent "dilution" of field enhancement. The 
growth of the positive streamer depends more non-linearly on the local field through 
impact ionization; therefore it stays thinner, the field is enhanced more and subsequently 
it propagates faster. In experiments, this asymmetry is found to decrease with increasing 
voltage; whether simulations show the same, will also have to be investigated in the 
future. 

Streamer physics is an exciting and widely open field, and it is amazing to note how 
little their polarity dependence has been characterized up to now, both experimentally 
and theoretically. Of course, the study of single streamers is only one problem in a range 
of phenomena, other questions concern streamer branching [42l |43] , interactions [24l |4T] , 
particle aspects [Ml US], the full inception process near an electrode or the electro- 
dynamical characterization of multi-streamer processes. 

Appendix A. Simulation of a needle electrode 

We simulated needle electrodes by a simplified version of the Charge Simulation 
Technique (CST) described in [16]. In general, the presence of an electrode imposes 
a fixed electrostatic potential along its surface. However, one can retain the main 
properties of a needle electrode by fixing the potential only at its tip (point P in Fig. lAll) . 
This is achieved by introducing a simulated point charge Q at a certain location inside 
the electrode. At each time step of the simulation, the value of Q is calculated to keep 
4>{P) (the electrostatic potential at P) fixed to 0(-P) = Vq. This schematic approach 
approximates the effect of a needle with a radius equal to the distance between P and 
Q {Rneedie in Fig- lAlj) and a length equal to the distance between the upper planar 
electrode and Q {Lneedie in Fig. jATj). 

In our simulations, we restrict the particles to the cylindrical volume below the 
needle tip and apply a homogeneous Neumann boundary condition at its top and bottom 
sides. Although this creates an artificial boundary inside the physical domain, note that 
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Figure Al. Implementation of a needle eleetrode using the Charge Simulation 
Technique (CST) with a single eharge. The particles are restricted to the rectangular 
domain below the needle, while the electrostatic field is solved in a larger domain 
limited by the two planar electrodes. 



our streamers will touch this plane only around the needle tip. Hence it can be used as 
a rough approximation for an electrode with a free in- or outflow of electrons. 



Appendix B. Electro-dynamic characterization of the streamer 

From one numerical simulation, one gets spatial profiles of the densities of species and 
electric fields. However, to help the physical understanding of the streamer process 
it is often useful to use more macroscopic quantities that roughly characterize the 
propagation of a streamer. Although the meaning of quantities such as radius or length 
is intuitively obvious, there is a certain unavoidable arbitrariness in the way one defines 
them from the more microscopic data. Therefore we provide here a precise description 
of the way we define each of these quantities. 

First of all we define the streamer tip as the point on the propagation axis where 
the absolute value of charge density is maximal. The space-charge layer of the streamer 
is then defined as the volume around the streamer tip where the absolute value of the 
charge density is larger than half of its maximum value. The streamer head charge Q 
is defined as the net charge content of the space-charge layer. 

The streamer length L is the separation between the needle electrode and the 
streamer tip. The streamer velocity is simply v = dL/dt. The enhanced field Emax is 
defined as the maximum of the electric field in absolute value. 

The definition of a streamer radius is somewhat more involved. We are mainly 
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interested in characterizing the shape of the space-charge layer. Hence we followed this 
procedure: for each z we took the radius r with the highest charge density. This gives 
us a r{z) curve that we restrict to the points inside the space-charge layer. This curve 
is fitted to a circle and we take the resulting radius as the radius of the streamer, R. 
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